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ABSTRACT 



Investigation of oceanographic data acquired in the Columbia River 
estuary, which lies between the states of Washington and Oregon, 
revealed significant variation in sound velocity with respect to time and 
location. National Ocean Service (NOS) standards for echo-sounder 
acquired soundings require the knowledge of sound velocity to within 
plus or minus 4 meters per second. To meet this requirement in the 
Columbia River estuary, methods were devised to allow prediction of 
sound velocity based on location and the height of tide. 

Confidence intervals associated with sound velocity predictions 
provided by regression analyses showed a substantial improvement in 
accuracy when compared to using a single average velocity for each 
location. Nevertheless, there were still some cases where the NOS 
requirement could not be met. 
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I. INTRODUCTION 



A. NOS HYDROGRAPHY 

"Hydrography is that branch of physical oceanography dealing with 
the measurement and definition of the configuration of the bottoms and 
adjacent land areas of oceans, lakes, rivers, harbors, and other water 
forms on earth" [Ref. 1: p. 1-3]. The National Ocean Service (NOS) 
is the federal agency responsible for producing accurate nautical charts 
of nearshore areas of the United States and its territorial possesions. 
The principal objective of most hydrographic surveys conducted by NOS 
is to obtain basic data for the compilation of nautical charts with 
emphasis on the features that may affect safe navigation [Ref. 1: p. 
1-3]. For this thesis data analysis is limited to the study of sound 
velocity variations and the correction of echo soundings for NOS 
hydrographic surveys. The type of data used for this investigation is 
often the basis for other physical studies such as the dynamics of 
oceanic mixing and salt transport. No attempt will be made to analyze 
or explain such processes. 

1 . Sound Velocity and Sound Speed 

Manuals and papers about hydrographic surveying often use 
the term "sound velocity" to designate what a physicist would define as 
"sound speed." The term "sound velocity" will be used in this thesis 
as many of the references also use the term. 

2. Echo Sounding 

Echo sounders are the basic tool of the NOS hydrographer for 
determination of depth. Echo sounders measure the time required for a 
sound wave to travel from its point of origin to the bottom and return. 
They also convert this measured time to distance or depth [Ref. 1: p. 
AG-1]. The basic formula for echo sounder depth determination is: 
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Z = VT/2 



( 1 . 1 ) 



where Z is the depth, V is the calibrated velocity of sound for the echo 
sounder, and T is the time between transmission and reception of the 
sound signal. For depths to be accurate precise measurements of time 
and a knowledge of the mean, or harmonic, sound velocity (Chapter II) 
are mandatory. 

B. THE VELOCITY OF SOUND IN SEA WATER 

Wilson [Ref. 2: p. 641] developed empirical equations based on 
laboratory measurements of sound velocity under tightly controlled 
conditions. The basic equation is: 

V = 1449.14 * Vj * Vp * Vg * VsTp (1.2) 

where V is the velocity of sound in meters per second, Vj is a 
correction for temperature, Vp is a correction for pressure, is a 
correction for salinity, and V^jp is a simultaneous correction for the 
combined effects of temperature, salinity, and depth [Ref. 1: p. 4-75]. 
Wilson's equations (Chapter II), which are currently used by NOS, 
were chosen to analyze the data for this thesis. 

C. NOS ACCURACY REQUIREMENTS 

The sound velocity at which NOS echo sounders are calibrated is 
800 fathoms per second. To obtain accurate depths for nautical charts 
the hydrographer must correct all soundings where the velocity of 
sound deviates from the calibrated value. There are several methods 
available for determining corrections for soundings, and they will be 
discussed in the following section. The NOS accuracy requirements 
[Ref. 1: p. 4-70] are: 

For use in correcting echo soundings, the velocity of sound must 
be known with sufficient accuracy to ensure that no sounding will 
be in error by as much as 0.25 percent of the depth from this 
cause alone. Therefore, the mean velocity must Be known to 
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within plus or minus 4 meters per second. To satisfy these 
accuracy requirements, one must know the mean temperature of the 
water to an accuracy of plus or minus 1 degree Celsius and the 
salinity to within plus or minus 1 part per thousand. A sufficient 
number of salinity and temperature observations must be made so 
that the velocity can be determined within the specified accuracy 
over the entire area sounded. 



D. PROCEDURES FOR CORRECTING SOUNDINGS 

Echo soundings can be corrected by either a direct comparison or 
an oceanographic procedure. Direct comparisons involve comparing echo 
soundings to known depths; this does not require a knowledge of the 
sound velocity. Procedures which involve the calculation of sound 
velocity from measurements of water characteristics are classified as 
"oceanographic and limnologic determinations" [Ref. 1: p. 4-74]. For 
this thesis the direct measurement of sound velocity by a sound 
velocimeter is also classified as an oceanographic procedure. 

1 . Direct Comparison Systems * 

Two procedures used to determine depth corrections by the 
direct comparison method are the bar check and the vertical cast. In 
these the effects of sound velocity variation and instrument error are 
included in the results [Ref. 1: p. 4-71]. 

a. Bar Check 

The bar check procedure involves lowering a bar over 
the side of a survey vessel and positioning the bar under the 
transducer at a known depth indicated by markings on the lines 
supporting the bar. Comparisons between the depth as indicated by 
the echo sounder and the known depth are made at 5- or 6-foot 
increments as the bar is lowered and raised though the water column. 
Due to the handling of the bar, this procedure is usually limited to 
smaller vessels and the quality of the results is a function of sea 
•conditions. Observations during rough seas are subject to unacceptable 
magnitudes of error [Ref. 1: p. 4-71]. 
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b. Vertical Cast 



The vertical cast requires simultaneous sounding of the 
bottom using an echo sounder and a lead line. This procedure, which 
can be accomplished by practically all survey vessels, is not as 
accurate as the bar check, but vertical casts can be used to determine 
or verify instrument errors. Where sounding sheets junction, there are 
often discrepancies and soundings fail to agree. Vertical casts can 
sometimes help resolve such problems [Ref. 1: p. 4-72]. 

2. Oceanographic Sound Velocity Determinations 

The more traditional oceanographic procedures involve sound 
velocity determination by measurement of the variables which enter 
Equation 1.2, namely temperature and salinity as functions of depth. 

The Nansen cast and electronic sensors such as the CTD (conductivity, 
temperature, and depth) or the STD (salinity, temperature, and depth) 
are examples of equipment used for these procedures. An alternative is 
the direct measurement of sdund velocity using a sound velocimeter. 

The electronic sensors record data almost continuously with depth. 

[Ref. 1: p. 4-74]. 

To obtain velocity corrections with Nansen, STD, or CTD 
data, Umbach [Ref. 1: p. 4-75] recommends using the "summation of 
layers" method, which extends from the concept of mean or harmonic 
sound velocity. 

Most of the electronic sensors which are currently in use yield 
precise measurements in excess of NOS accuracy requirements for sound 
velocity corrections if properly calibrated and maintained [Ref. 1: p. 
4-74] . According to NOS standards all electronic sensors must be 
periodically checked for accuracy against Nansen cast or equivalent 
observations . 

E. SOUND VELOCITY VARIATION 

The effect of changes in temperature, salinity, and pressure on 
sound velocity depends upon the initial conditions of these variables. 
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Wilson’s equations are used to illustrate in Table I the effect of 
changing variables by one unit. 



TABLE I 

Variation of Sound Velocity 



PRES 

db 


TEMP 


SAL 

g/kg 


V 

m/s 


delta 

m/s 


1.00 

1.00 


15.00 

14.00 


30.00 

30.00 


1502.90 

1499.30 


3.60 


1.00 

1.00 


5.00 

4.00 


30.00 

30.00 


1464.27 

1460.03 


4.25 


1.00 

1.00 


10.00 

10.00 


35.00 

34.00 


1490.76 

1489.48 


1.28 


1.00 

1.00 


10.00 

10.00 


1.00 

0.00 


1449.01 

1447.84 


1.17 


60.00 

59.00 


10.00 

10.00 


30.00 

30.00 


1485.33 

1485.32 


0.02 


1.00 
0.00 . 


10.00 

10.00 


30.00 

30.00 


1484.38 

1484.36 


0.02 



Velocity (V) was computed using Wilson's equations. 



Corrections for .sound velocity variations in estuaries can present 
the hydrographer with a difficult task. According to Cameron and 
Pritchard [Ref. 3: p. 306], an estuary is a semi-enclosed coastal body 
of water having a free connection with the open sea and within which 
the sea water is measurably diluted with fresh water deriving from land 
drainage. Estuaries are extremely dynamic systems. Topography, 
river flow, and tidal action are important factors that influence the rate 
and extent of mixing salt and fresh water [Ref. 4: p. 4]. For the 
hydrographer, estuaries are areas where sound velocity varies not only 
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with location but also in time due to the tidal cycle and variations in 
river flow. Although no specific procedures are given, Umbach 
[Ref. 1: p. 4-75] notes the possible complexity of such areas: 



Under special circumstances such as ... in areas of extensive 
estuarine discharge of fresh water . . . series of regional velocity 
curves are required. Regional curves must be carefully studied to 
determine how they can be best grouped by area or time or to 
permit drawing of average curves. 



F. OBJECTIVES 

One objective of this thesis was to investigate the magnitude of 
sound velocity variation associated with significant changes in water 
characteristics within an estuary. 

A second objective was to investigate the feasibility of determining 
velocity corrections according to NOS standards in areas where extreme 
changes in water characteristics occur with respect, to space and time. 

A third objective was to write a FORTRAN program to enable 
further testing by the field hydrographer equipped with a personal 
computer. The use of commercial subroutines has been avoided. 
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II. METHODOLOGY AND PROCEDURE 



A. AREA OF STUDY 

The Columbia River estuary which discharges into the Pacific 
Ocean at latitude 46°15'N, longitude 124°00'W has variations in water 
characteristics which make it an excellent area for investigation. 
According to Dyer [Ref. 4: p. 96], the Columbia River estuary is 
characterized by an unusual combination of large river discharge and 
strong tidal currents. Recent investigations have revealed differences 
in salinity of up to 30 grams per kilogram at individual stations during 
a single tidal cycle [Ref. 5: p, 104]. Due to this large difference and 
the availability of data from a recently completed survey, this region 
was chosen for investigation. 

B. DATA SOURCES 

To analyze the problem completely, water characteristics 
(temperature and salinity as functions of depth) must be known at 
stations spaced throughout the region, and each station occupied during 
a complete lunar tidal cycle. Variations caused by seasonal changes in 
fresh-water runoff should be examined thoroughly. 

1 . Columbia River Estuary Data Development Program 

The Columbia River Estuary Data Development Program 
(CREDDP) was a 5-year, federally funded study which included 
approximately 2 years of data acquisition for investigation of physical 
and biological characteristics. CREDDP was designed to meet the needs 
of two major groups. One group consisted of local, state, and federal 
government agencies involved in planning and permitting activities. 

The second group included research scientists and educational 
institutions [Ref. 6: p. v] . 
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One accomplishment of CREDDP was the acquisition of current 
velocity, conductivity, temperature, and depth (VCTD) data for a 
better understanding of physical processes such as currents, sediment 
transport, and salinity intrusion. VCTD profiling was carried out for 
13 consecutive days from 17 stations for periods of from 12 to 48 hours 
per station. Between October 16 and 27, 1980, approximately 1600 
VCTD casts were made. Of these, 458 casts were made available to 
users [Ref. 6: p. L-8] . 

Data were supplied to NPS by the U.S. Army Corps of 
Engineers (COE) at Portland, Oregon, on two 9-track magnetic tapes 
[Ref. 6: pp. 1-5]. 

For an analysis of sound velocity variation, VCTD casts from 
11 of the 17 stations were investigated. The stations, along with their 
positions and times are listed in Table II. The plot of the geographical 
positions of the stations is shown in Figure 2.1. These stations 
covered the area where the maximum sound velocity variation would be 
expected as deduced from salinity intrusion diagrams of the estuary 
[Ref. 5: p. 104]. 

2. NOS Bar Check Data 

NOS bar check data were used for a sound velocity variation 
analysis. The descriptive report for survey H-8422 listed the day each 
bar check was taken and the correction for each 5-foot increment of 
depth. Also given was a general term for the stage of the tide such as 
ebb, flood, high water, or low water. Data for 46 bar checks were 
acquired between May 9 and July 16, 1958 [Ref. 7: p. 1]. 

3. NOS Predicted Tide Data 

NOS predicted tide data were obtained for correlation and 
regression with the CREDDP and the H-8422 bar check data. Predicted 
tides were used because field hydrographers rarely have the time or 
the proper tools for a detailed analysis of real tide data. 
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table II 

Columbia River Estuary Data Development Program 



STATION (STN) INFORMATION 



STN 


POSITION 


ACOUISITION 

START 


TIMES (PDT) 
END 




2N 


46°15' 
124°03 ' 


'45"N 

'31"W 


10/18/80 


0930 


10/19/80 


0230 


2S 


46°15' 
124°03 ' 


' 13"N 
' 20"W 


10/18/80 


1230 


10/19/80 


0930 


4NA 


46°14' 

123°56' 


03"N 

'09"W 


10/20/80 


2115 


10/21/80 


1745 


4NB 


46°14' 

123°55' 


' 47 "N 
' 53"W 


10/20/80 


2105 


10/21/80 


1730 


4SA 


46°12' 

123°57' 


'42 "N 
09"W 


10/26/80 


0007 


10/27/80 


0832 


4SB 


46°12' 

123°56’ 


' 51"N 
'55"W 


10/25/80 


2350 


10/27/80 


0820 


SNA 


46°13' 
123°53 ' 


'55"N 
' 17 "W 


10/20/80 


1733 


10/22/80 


0930 


5NB 


46°14' 

123°53' 


' 18"N 
’ 14"W 


10/20/80 


2035 


10/22/80 


0913 


5NC 


46°14'07"N 
123°53 ' 14"W 


10/21/80 


0935 


10/22/80 


1003 


5SA 


46°11'55"N 

123°48'38"W 


10/24/80 


0830 


10/25/80 


2331 


5SB 


46°12 '05"N 
123°48 ' 54"W 


10/24/80 


1126 


10/25/80 


1130 
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COLUMBIA RIVER ESTUARY 

CREDDP STATIONS 




Figure 2.1 CREDDP Stations. 

C. CREDDP DATA ANALYSIS 
1. COE V CT D Data Tapes 

COE VCTD data tapes contained pertinent information such as 
header data, time, temperature, depth, conductivity, and salinity. 
Mercier [Ref. 6: p. L-8] shows a breakdown of the data components for 
one record of data. The tapes contained more data records than were 
necessary for this investigation. The VCTD sampling device recorded 
data every 0.04 seconds, thereby making the average number of 
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records for a single cast approximately 1500. Having 270 casts to 
analyze, a method was devised to reduce the data to a more efficient 
size. A computer program was formulated to input the data for each 
cast and then output the first record and the records which 
corresponded to the nearest whole meter of depth. For example, a 
10-meter cast would theoretically have 11 records of data; the first 
record would represent the surface layer, and the remaining 10 records 
would be for 1 -meter increments. In accordance with the NOS policy of 
using bar check data measured at 5- or 6-foot increments, data were 
eliminated except at 1 -meter increments. A sample of a profile after 
editing is shown in Table III. Standard depths which represent depths 
of 1, 3, 5, 7, 9, 11, and 13 meters, were adopted to further simplify 
data processing and description. 

Several drawbacks were encountered using the CREDDP data. 
Because of the shallow depths of CREDDP stations, many of the casts 
were inadequate for correcting soundings in the deeper areas of the 
estuary. Some casts did not include surface layer data (0 to 2 
meters) . 

2. Sound Velocity Computations 
a. Wilson's Equations 

Sound velocity was computed for each record by Wilson's 
equations which use temperature, pressure, and salinity for the 
computation. In the following equations V is sound velocity in meters 
per second; T is temperature in degrees Celsius, where -4 < T <30 
°C; P is pressure in kilograms per square centimeter, where 1 < P < 
1000 kg/cm'^' and S is salinity in grams per kilogram, where 0 < S < 37 
g/kg. 

V = 1449.14 + Vj + Vp Vs * Vsjp (2.1) 

Vj = 4.5721T - 4.4532x10’2t2 - (2.2) 

2. 6045x1 O'^T^ * 7.9851 xIO'^T^ 
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Vp = 1. 60272x10" V * 1 .0268xl0"^p2 + (2.3) 

3.5216xl0"^p2 - 3.3603x10"''2p4 

Vg = 1.39799(3-35) ♦ 1 .69202xl0"^(S-35)2 (2.4) 

Vg-pp = (S-35)(-1.1244 x10"2t * 7.7711x10"^t2 + (2.5) 

7.7016xl0"^P - 1.2943xl0"^p2 ♦ 

3.1580xl0"^PT * 1.5790 x10"^Pt2) * 

P(-1.8607xl0"'^T * 7.4812 x10"®t2 * 

4. 5283x1 0"^T^) * p2(-2.5294xl0"^T + 

1.8563x10"^t 2) * p3(-l. 9646x10" ''^T) 

The standard deviation for the velocity is 0.30 meters per 
second [Ref. 8: p. 1357]. 







TABLE 


III 








4SA 


--Cast 6 


--Sound Velocity 


Variables 








COLUMBIA RIVER 








STATION = 


4SA„ , 




CAST = 031B08 




LATITUDE 


46®12' 


42 "N 


LONGITUDE 123®57'09 


"W 


DATE OF CAST 1980 


OCT 26 


TIME 


; OF CASl 


' 5:16: 


0 






CTD 










TIME 


PRES 


DEPTH 


TEMP 


COND 


SAL 




s 


db 


m 


°C 


ratio 


g/kg 




0.04 


-0.04 


0.10 


12.26 


0.4531 


15.71 




5.48 


0.85 


1.00 


12.04 


0.5252 


18.57 




12.04 


1.86 


2 . 01 


11.79 


0.7073 


25.91 




18.48 


2.85 


3.00 


11.78 


0.7571 


27.95 




24.60 


3.86 


4.01 


11.78 


0.7932 


29.42 




31.52 


4.86 


5.00 


11.79 


0.8184 


30.46 




37.68 


5.90 


6.04 


11.79 


0.8164 


30.37 




43.80 


6.86 


7.00 


11.80 


0.8232 


30.64 




50.48 


7.88 


8.01 


11.81 


0.8291 


30.88 




56.80 


8.89 


9.01 


11.81 


0.8302 


30.92 




63.32 


9.88 


10.00 


11.81 


0.8330 


31.04 




69.92 


10.91 


11.03 


11.82 


0.8338 


31.07 




76.12 


11.90 


12.01 


11.82 


0.8356 


31.14 




82.84 


12.90 


13.01 


11.82 


0.8346 


31.10 
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b. Mean Sound Velocity 

Because echo sounders are time measuring devices, 
distance or depth is one-half the round trip travel time of a sound 
pulse multiplied by a suitable mean velocity. The mean sound velocity 
(MSV) is given by 

MSV = Z/T (2.6) 

where Z is the distance traveled in time T. To correct for sound 
velocity variations, the water column is divided into a series of layers, 
each with thickness AZj and an associated velocity Vj. Therefore, the 
time interval for the sound wave to pass through the i^*^ layer is: 



ATj = AZj/Vi (2.7) 

Summing all layers results in; 

T = 2 ATj =2 AZ/V: (2.8) 

i=l ' i=l ' ' • 

or; 

T = /odZ/V (2.9) 

Substituting into Equation 2.6 results in: 

MSV = Z/(/^dZ/V) (2.10) 



Equation 2.10 was evaluated using the trapezoidal method [Ref. 9: p. 
86]. NOS uses mean sound velocities in the "summation of layers" 
method described by Umbach [Ref. 1: p. 4-74] and Wallace [Ref. 10: 
p. RK530.1]. 

Two sound velocity listings computed using the data in 
Table III are shown in Table IV. The sound velocity column was 
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computed using Equation 2.1, and the mean velocity column was 
computed using Equation 2.10. 



TABLE IV 

4SA--CAST 6--Sound Velocity Data 



DATE: 


26 OCTOBER 1980 


TIME: 


0520 


(PDT) 




SOUND 


MEAN 


DEPTH 


VELOCITY 


VELOCITY 


m 


m/s 


m/s 


0.10 


1475.5 


1475.5 


1.00 


1478.1 


1476.8 


2.01 


1486.1 


1479.6 


3.00 


1488.6 


1482.2 


4.01 


1490.5 


1484.1 


5.00 


1491.8 


1485.5 


6.04 


1491.7 


1486.6 


7.00 


1492.1 


1487.3 


8.01 


1492 . 5 


1488.0 


9.01 


1492.5 


1488.4 


10.00 


1492.7 


1488.9 


11.03 


1492.8 


1489.3 


12.01 


1492.9 


1489.5 


13 . 01 


1492.9 


1489.8 


14.01 


1492.9 


1490.0 



To simplify the analysis of the variation of mean sound 
velocity as a function of time, a computer program was written to 
arrange the data in individual arrays for each depth at each station. 

3 . Mean Sound Velocity and Tide Height 

Inspection of the mean sound velocity arrays suggested that 
velocities were maximum at high tide and minimum at low tide. To test 
this, a correlation and regression analysis between the CREDDP derived 
mean sound velocities and predicted tide heights was performed. 
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a. Predicted Tide Heights 

To perform the correlation and regression analysis, a 
method had to be developed to calculate the height of the tide at any 
desired time. One method would be to obtain the harmonic constants 
for prediction of the Columbia River estuary tides from the NOS Tides 
and Water Levels Branch. These constants would then be used in the 
appropriate equations [Ref. 11: p. 123]. From a computer programming 
standpoint this method would be the simplest. A second method would 
be to formulate a computer program based on the graphical method 
(Figure 2.2) described in the NOS Hydrographic Manual [Ref. 1: p. 
4-68]. This method was chosen because normally the field 
hydrographer is not supplied with the harmonic constants for project 
areas. However, the tide tables are an essential tool of the 
hydrographer . 

The present NOS HYDROPLOT computer program for tide 
prediction- lacks documentation and is not designed to solve this 
particular problem [Ref. 10: p. 500.1]. Therefore, a program was 
written. 

Points A and E (Figure 2.2) represent the predicted low 
and high tide, respectively. The straight line between these points 
was divided into four equal parts by points B, C, and D. Point B' was 
plotted directly bejow B at one-tenth the total range, and point D' was 
plotted directly above D using the same distance [Ref. 1: p. 4-68], 
forming three regions for the solution of predicted heights between 
points A and E. For the region between points A and B’, the curve 
was assumed to be parabolic, and in order to calculate the tide height 
for a designated time in this region, the equation 

Hj = aTj2 + bTj + c (2.11) 

was used, where H is the tide height and T is the designated time. To 
calculate the coefficients a, b, and c; the equations 

Ha = ^ (2.12) 
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Figure 2.2 Predicted Tide Graph. 
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c 



(2.13) 



Hgt = aTg'^ ♦ * 

0 = 2aT^ ♦ b (2.14) 

were used. For times that lie between points D' and E, the same 
method was used. Between points B' and D' the change in height was 
considered to be linear with respect to time and can be computed by 
the equation 

Hj = Hg. * m(Tj-Tg-) (2.15) 

where m is the slope defined by 

m = (Hp-HgO/(Tp.-Tg.) (2.16) 

Using the program based on the above equations, tide heights were 
calculated for the times of all VCTD casts used in the analysis. "Data 
Array" means a set of sound velocity and predicted tide height data 
that shows the variation of both variables with time for a particular 
station and depth (Table V). 

b. Data Plots 

Computer generated plots of predicted tide height and 
observed mean sound velocity (Figure 2.3) were made for all stations 
using mean velocities based on standard depths. 

c. Correlation Analysis 

A correlation analysis routine for mean sound velocity 
and tide height was run for all data arrays. In addition to computing 
the correlation coefficient, the routine also computed the averages and 
the standard deviations of the mean sound velocity and the predicted 
tide height for each data array. Averages were computed by the 
equation 

7 = :S,X;/n (2.17) 

1 = 1 ' 
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4SB--Data 


TABLE 

Array-- 


V 

•Depth = 7 m 




TIDE 


MEAN 


TIME 


HEIGHT 


VELOCITY 


hr 


ft 


m/s 


0.12 


1.2 


1469.6 


1.50 


4.1 


1476.6 


2.63 


6.5 


1486.5 


3.60 


7.6 


1490.8 


4.62 


7.6 


1490.8 


5.27 


7.1 


1487.3 


6.32 


5.6 


1483 . 1 


9.23 


1.7 


1472 . 1 


11.55 


3.0 


1471.0 


12.58 


4.9 


1472.2 


13 . 17 


6.0 


1476.9 


17.53 


7.7 


1489.6 


21.52 


-0.8 


1471.4 


22.43 


-1.5 


1468.6 


23.35 


-1.3 


1467.7 


24.28 


-0.2 


1467 . 1 


25.33 


1.9 


1470.0 


26.42 


4.1 


1473 . 6 


27.57 


6.2 


1480.0 


28.48 


7.2 


1485.9 


29.48 


7.2 


1491.9 


30.50 


6.4 


1488.6 


31.53 


5.1 


1485.4 


32.53 


3.8 


1477.3 



where Xj is the observation and n is the number of observations 
[Ref. 12: p. 124]. Standard deviations [Ref. 12: p. 128] were 
computed using the formula 

s = /s (Xj-x)^/(n-1 ) (2.18) 

V 1 = 1 

Computation of the simple correlation coefficient [Ref. 12: p. 258] of 
the mean sound velocity and predicted tide height arrays was done with 
equation 



r = [n XjYj - ( J^j)(,|:T Yj)] 

/ n ^ n rx n o 

/[n2^Xi2 - (.2 Yi^)) 



(2.19) 
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Figure 2.3 4SB--Data Array Plot--Depth = 7 m. 
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A sample of the results of a correlation analysis is listed in Table VI. 



TABLE VI 






4SB--Statistical Data--Depth = 7 m 






AVERAGE TIDE HEIGHT 


4.4 


ft 


STANDARD DEVIATION OF TIDE HEIGHT 


3.0 


ft 


AVERAGE MEAN VELOCITY 


1480.8 


m/s 


STANDARD DEVIATION OF MEAN VELOCITY 


8.0 


m/s 


CORRELATION COEFFICIENT 


0.93 



d. Regression Analysis 

To determine the feasibility of predicting mean sound 
velocity, a first-order (or linear) regression analysis was performed on 
all data arrays. The data were fitted to the equation 

y = a + bx (2.20) 

where y is the estimated mean velocity, a is the coefficient, b is the 
the intercept, and x is the tide height. The coefficient and intercept 
estimates were determined to minimize the sum of the residuals 

|l®i = litYi-Vi) (2.21) 

where Yj is the observed mean sound velocity and yj is the estimated 
velocity [Ref. 13: p. 533]. Determination of the coefficients and the 
intercepts was accomplished by the normal equations 



28 



( 2 . 22 ) 



n ri 

an + b 2 X: = 2 Y: 

1=1 ' i=l ' 

n n o n 

a 2 X: * b.2 Xj2 = 2 X:Y: (2.23) 

i=l ' i=l ' i=l ' ' 

where n is the number of observations for a particular station and 
depth [Ref. 13: p. 533]. The sum of the squares of the residuals is 

SS = 2 (Yj-Yj)2 (2.24) 

The standard error of estimate is 

see ^(SS / v) (2.25) 

where v is the degrees of freedom given 

V = p - n (2.26) 

where n is the number of sample observations and p is the number of 
estimated variables [Ref; 14: p. 131]. For a first-order analysis, there 
are two estimated variables. Results of a first-order analysis are listed 
in Table VII. 

Confidence intervals associated with predicted mean 
velocities obtained by equations derived from the linear regression 
analysis were calculated to determine the. accuracy of predictions. 
Confidence interval computation requires the standard error of estimate, 
the value of the student-t (t) distribution which is based on the 
desired percent of confidence, and the degrees of freedom for the 
sample. For this thesis, 90-percent confidence intervals were computed 
using the equation 

Cl = ± t(v,0.95)[1 Ml/n) * (2.27) 

n 

[(Xj-x)^ / 2 (xj-x)^]]see 
i = l ' 
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TABLE VII 

4SB--Linear Regression--Depth = 7 m 



COEFFICIENT = 
Y INTERCEPT = 



CAST 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

SS 

see 



1468 !0 



MEAN 

VELOCITY 

RESIDUAL 



- 0 , 
0 , 
0 , 
3. 
2 . 
0 . 
-5 



- 0.8 

-1.7 



' 1 , 

■ 1 , 

- 1 , 

0 . 

■ 0 , 



2 

0 

4 

1 

9 



-4.5 

-4.0 



1 

5 

7, 

5 

2 



4 
0 
2 

5 
9 



0.0 

-3.2 

- 1.2 

0.0 

5.1 

4.8 

0.0 

-1.7 

- 0.8 

244.8 

3.0 



A RESIDUAL IS THE DIFFERENCE BETWEEN THE OBSERVED 
MEAN VELOCITY AND THE PREDICTED MEAN VELOCITY. 



where Cl is the confidence interval for a single estimated velocity (yj), 
Xj is the tide height, "x is the mean tide height, and n is the number of 
sample observations [Ref. 15: p. 30]. 
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Several plots of the data arrays indicated a curvilinear 



relationship between mean sound velocity and predicted tide height. 
Therefore, second- and third-order polynomial regression analyses were 
performed on all data arrays to provide a better fit of the data, 
thereby increasing the accuracy of the mean sound velocity prediction. 
The normal equation method used in the first-order analysis was 
extended for calculation of a second-order polynomial. To check the 
results of the programs for the normal equation method, a subroutine 
using the orthogonal polynomial method was obtained from the NPS 
Computer System Library [Ref. 16: p. 8 ]. A third-order regresssion 
analysis was also performed by the orthogonal polynomial subroutine to 
investigate the possibility of increasing the accuracy of prediction and 
fit. Based on results of the investigation and check (Chapter III), the 

A 

first-order and second-order normal equation methods were designated 
as the primary methods of analysis. The normal equation method 
appeared to be the simplest to program. The second-order analysis 
uses the equation 

y = ag aj^ X + a 2 (2.28) 

where y represents the estimated mean sound velocity and x is the 
predicted tide height [Ref. 13: p. 535]. Coefficients (ag, a-^ , and 82 ) 
were determined to minimize the sum of residuals by the following 
normal equations 



where n is the number of casts for a particular station and depth 
[Ref. 13: p. 536]. Results of a second-order analysis are listed in 
Table VIII. 




(2.29) 




(2.30) 




(2.31) 
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TABLE VIII 

4SB--Polynotnial Regression--Depth = 7 m 



COEFFICIENTS 





an 


= 1470.0 




aV 


= 1.6 




a2 


= 0.1 






MEAN 

VELOCITY 


CAST 




RESIDUAL 


1 




-1.7 


2 




0.8 


3 




1.0 


4 




3.2 


5 




2.7 


6 




0.2 


7 




-4.a7 


8 




0.6 


9 




-1.0 


10 




-1.7 


11 




-1.4 


12 




-0.9 


13 




-0.7 


14 




-2.4 


15 




-5.0 


16 




-3.2 


17. 




2.0 


18 




4.2 


19 




5.3 


20 




1.4 


21 




-0.6 


22 




-1.7 


23 




-3.0 


24 




-0.4 


25 




0.8 


26 




5.4 


27 




5.0 


28 




0.7 


29 




-0.9 


30 




0.0 


SS 




209.8 


see 




2.8 



A RESIDUAL IS THE DIFFERENCE BETWEEN THE OBSERVED 
MEAN VELOCITY AND THE PREDICTED MEAN VELOCITY. 



Confidence intervals for predicted velocities were 
computed using Equation 2.27. The degrees of freedom for a 
second-order analysis was less than a first-order analysis because of 
the addition of one estimated variable. 
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D. NOS BAR CHECK DATA 



The NOS bar check data consisted of 46 bar checks which were 
acquired to determine the corrections for soundings of survey H-8422. 
The format of the data was a listing of corrections for 5-foot increments 
of depth. The initial analysis included determination of the maximum, 
minimum, and mean correction for each depth (Table IX). 



Bar 


Check 


TABLE IX 

Correction Data For Survey H 


-8422 


OBS 


MEAN 

ft 


UPPER 

LIMIT 

ft 


LOWER 

LIMIT 

ft 


RANGE 

ft 


DEPTH 

ft 


88 


-1.0 


-0.6 


-1.4 


0.8 


5 


90 


-0.8 


-0.2 


-1.4 


1.2 


10 


90 


-0.7 


-0.1 


-1.4 


1.3 


15 


90 


-0 . 7 


.0.0 


-1.3 


1.3 


20 


90 


-0.7 


0. 1 


-1.4 


1.5 


25 


87 


-0.7 


0.0 


-1.6 


1.6 


30 


. 82 


-0.6 


0.0 


-1.3 


1.3 


35 


73 


-0.5 


0.0 


-1.2 


1.2 


40 


30 


-0.3 


0.6 


-1.1 


1.7 


45 


12 


-0.4 


0.0 


-1.0 


1.0 


50 



For each depth of a bar check there were two corrections listed; 
the first correction was recorded as the bar was being lowered, and the 
second during raising. These two corrections were averaged to obtain 
a single value which was used for a correlation analysis with predicted 
tide height. 
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A plot of the bar check and tide height data (Figure 2.4) was 
made and visually examined. The plot indicates a poor correlation 
between bar check corrections and predicted tide heights. For 
verification, a correlation analysis (Table X) was performed using 
methods previously discussed. 



TABLE X 

Bar Check Statistical Data for Survey 


H-8422 




MEAN TIDE HEIGHT 


. . 3.3 


ft 


STANDARD DEVIATION OF TIDE HEIGHT 


. . 2.9 


ft 


MEAN CORRECTION 


.. -0.7 


ft 


STANDARD DEVIATION OF CORRECTION 


. . 0.3 


ft 


CORRELATION COEFFICIENT 


.. 0.12 




DEPTH OF OBSERVATIONS 


. . 25.0 


ft 



Because of the poor correlation, actual tidal data were obtained to 
investigate the possiblilty that the results were related to the predicted 
tides. No improvement in the results was accomplished by using actual 
tide data. No further manipulation of the NOS bar check data was 
carried out because of the number of unknowns associated with the 
acquisition procedures and because the amount of speculation necessary 
could not lead to meaningful analysis. 
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Figure 2.4 Bar Check Data Array--Depth = 25 ft. 
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III. DISCUSSION OF RESULTS 



A. CREDDP DATA 

Data acquisition procedures of CREDDP made an investigation of 
mean sound velocity variation possible, but there are factors which 
make the study incomplete. Because of this, a complete model for mean 
sound velocity prediction and correction based on methods presented in 
this thesis was not fully developed. However, the primary components 
of such a model were devised and will be explained in this chapter. 

The factors which make the investigation incomplete and 
recommendations concerning verification of ideas and methods presented 

A 

in this thesis will be addressed in Chapter IV. 

Investigation of the CREDDP data revealed that sound velocity in 
the Columbia River estuary varies with position, depth, and time. To 
analyze the problem the study area was divided into two primary 
sections, the north channel and the south channel (Table XI). The 
'river mile' distance represents the shortest deep-water path from buoy 
"7" at the river entrance to the station. The south channel, which is 
maintained for commercial marine traffic, has a river mile scale which 
can be found on NOS chart 18521, whereas north channel stations were 
scaled using chart 18521. 

Two categories of mean sound velocity variation are considered. 
The term "horizontal variation" designates differences in mean sound 
velocity between stations; i.e., how the velocity changes with location. 
The horizontal variation analysis is needed to develop zones or regions 
where velocity correctors may be grouped. The term "temporal 
variation" designates the change in mean sound velocity with time at 
each station. With temporal variation, mean sound velocity correctors 
for each zone change with time. 
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TABLE 

Channel 


XI 

Data 




north 


channel 


south 


channel 


Station 


River Mile 


Station 


River Mile 


2N 


2.0 


2S 


2.0 


4NB 


7.5 


4SB 


8.0 


5NB 


9.5 


5SA 


15.5 



1 . Horizontal Variation 

For the analysis of horizontal variation the two requirements 
made were: (1) that for each depth there be but one horizontal 
variation analysis and; (2) the average of all mean velocities obtained at 
the depth being analyzed be designated as the average velocity. Use 
of the average velocity allows each station to have a single velocity for 
each depth. The average velocity presents a problem if more 
observations were made on one end of the tidal spectrum than the 
other. For example, if more casts were acquired during periods of 
high water than periods of low water, the average velocity would be 
higher than if the majority of casts were acquired during low water 
periods. One of the CREDDP objectives was to acquire casts on an 
hourly basis. This objective was met in most cases, but there were 
some instances where the time interval between successive casts was 
greater than 1 hour. For example, station 4SB was chosen over station 
4SA for designation of the south channel because there were no casts 
acquired at station 4SA during the higher high tide. Another potential 
problem would be introduced if there were significant variations in the 
river flow. According to Jay [Ref. 5: p. 114], river flow remained low 
for the acquisition period. An example of horizontal variation of the 
average velocity is shown in Figure 3.1. These plots indicate the rate 
of horizontal variation of the average velocity is approximately the same 
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for the north and south channels. Data for the 11- and 13-meter 
depths were limited and did not cover the range of tide well in many 
cases . 

Regional zoning of mean sound velocity, correctors based solely 
on horizontal variation can be accomplished with little difficulty. Using 
the depth with maximum horizontal variation, zones can be drawn based 
on the NOS criteria of plus or minus 4 meters per second. A major 
problem still exists in that the temporal variation of mean sound velocity 
has not been accounted for when only the horizontal variation of the 
average velocity is considered. 

2. Temporal Variation 

An analysis of temporal variation of the mean sound velocity 
demonstrated that proper correction “'of soundings in the Columbia River 
is a multi-dimensional problem. Sound velocity variables change 
significantly with time at the same location, and the magnitude of 
variation can depend on a number of cyclic and noncyclic phenomena. 
Sound velocity variables change with the tide, and the magnitude of the 
variation is affected by monthly tidal, events such as spring and neap 
tides. The temporal and horizontal variations of velocity for the 
stations of the north and south channels are shown in Figure 3.1. 

Data for stations 2N and 2S were acquired during a neap tide 
period when the predicted tide range was minimum, whereas data 
acquisition for stations 4SA, 4SB, 5SA,and 5SB was accomplished 
during a spring tide period. Data for stations 4NA, 4NB, 5NA, 5NB, 
and 5NC were acquired during a period of average tide range. As 
illustrated by Figure 3.1, stations 4NB and 4SB have about the same 
average velocity at the standard depths, but with respect to temporal 
variation, the range of velocity for station 4SB was significantly larger 
than that found at station 4NB. According to Jay [Ref. 5: pp 
114,117], the neap tide period was a one of minimum vertical mixing and 
maximum stratification, whereas during the spring tide period, maximum 
mixing, minimum stratification, and maximum salinity intrusion occurred. 
Examination of the casts for stations 4NB and 4SB (Table XII) 
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Figure 3.1 Temporal and Horizontal Variations of Velocity-- 

Depth = 9 m. 
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confirmed the previous statement. The cast for station 4NB acquired at 
predicted low tide showed a significant difference in temperature when 
compared to the low tide cast of station 4SB. The warmer temperatures 
of the low tide cast caused the the minimum velocity to be higher at 
station 4NB than at station 4SB. Comparison of the high tide casts 
indicated the surface layer salinity was 6.82 parts per thousand less at 
station 4NB than at station 4SB. The relatively large salinity 
difference for the high tide casts caused the maximum mean velocity to 
be lower at station 4NB when compared to the maximum mean velocity at 
station 4SB. Therefore, the increased energy associated with a larger 
tide range can increase the mean sound velocity range for certain areas 
in two ways. First, increased mixing associated with the larger tide 
range can moderate temperature for low tide cast, thereby causing 
minimum mean velocities to be lower. Second, increased salinity 
intrusion associated with a larger tide range can cause maximum mean 
velocities to be higher when compared to maximum mean velocities 
observed during an average or neap tide range. 

a. Correlation of Tide Height and Sound Velocity 

The effort to correct mean sound velocity for temporal 
variation began with a correlation analysis of predicted tide height and 
mean sound velocity. The results (Table XIII) indicate a strong 
relationship between tide height and mean sound velocity for most of 
the CREDDP data. Correlation in the 1-meter layer was usually not as 
strong as the deeper layers for most of the stations. Dashed line 
entries in Table XIII indicate insufficient data for a correlation 
analysis . 

b. Regression Analysis 

First- (linear), second-, and third-order polynomial 
regression analyses were performed for each standard depth of the 
CREDDP sampling stations. The objective of these analyses was to 
provide a numerical method of predicting mean sound velocity using the 
predicted tide height. These methods would be considered completely 
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TABLE XII 

The Effect of Tide Range on Mean Sound Velocity Range 



STATION 4NB 

TIDE HEIGHT = 0.02 ft (LOW TIDE) 



DEPTH 


TEMP 


SAL 


VEL 


m 


°C 


g/kg 


m/s 


1.0 


14.22 


6.32 


1472.2 


3.0 


14.14 


7.48 


1472.5 


5.0 


14.05 


8.52 


1473 . 1 


7.0 


13.93 


9.73 


1473.5 


TIDE HEIGHT = 8. 


3 ft (HIGH TIDE) 


DEPTH 


TEMP 


SAL 


VEL 


m 


°C 


g/kg 


m/s 


1.0 


12.18 


23.49 


1484.6 


3.0 


11.78 


29.26 


1487.4 


5.0 


11.97 


29.89 


1489.2 


7.0 


11.66 


30.19 


1489.8 




STATION 4SB 




TIDE HEIGHT = -1 


.5- ft 


(LOW TIDE) 


DEPTH 


TEMP 


SAL 


VEL 


m 


°C 


g/kg 


m/s 


1.0 


12.71 


8.41 


1468.6 


3.0 


12.69 


8.92 


1468.9 


5.0 


12 . 66 


9.55 


1469.2 


7.0 


12 . 64 


9.70 


1469.4 


TIDE HEIGHT = 9. 


2 ft (HIGH TIDE) 


DEPTH 


TEMP 


SAL 


VEL 


m 


°C 


g/kg 


m/s 


1.3 


11.88 


30.31 


1491.9 


3.0 


11.87 


30.69 


1492 . 1 


5.0 


11.84 


31.35 


1492 . 5 


7.0 


11.84 


31.48 


1492 . 7 



successful if the confidence interval for predictions based on all 
possible tide heights was less than, or equal to, plus or minus 4 meters 
per second. The employment of three different orders of polynomial 
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TABLE XIII 

Correlation Coefficients 



STN 


Im 


3m 


5m 


DEPTH 

7m 


9m 


11m 


13m 


2N 


0.61 


0.81 


0.85 


0.87 


0.87 


0.87 


0.78 


2S 


0.86 


0.68 


0.79 


0.83 


0.88 


0.89 


0.85 


4NA 


0.83 


0.91 


0.94 


0.97 


0.98 


0.99 


0.99 


4NB 


0.81 


0.90 


0.94 


0.96 


0.96 


0.99 


— 


4SA 


0.71 


0.82 


0.86 


0.89 


0.92 


0.90 


0.93 


4SB 


0.75 


0.83 


0.89 


0.93 


0.96 


0.99 


— 


SNA 


0.86 


0.89 


0.92 


0.91 


0.90 


— 


— 


5NB 


0.80 


0.86 


0.91 


0.94 


0.95 


0.97 


0.96 


5NC 


0.90 


0.91 


0.91 


0.94 


0.96 


0.97 


0.97 


5SA 


0.89 


0.88 


0.88 


0.87 


0.85 


0.86 


— 


5SB 


0.88 


0.88 


0.88 


0.87 


0.86 


— 


— 



regression analysis was for comparison of fit, as well as comparison of 
the predicted velocities and associated confidence intervals. 

A potential source of error may exist when using 
predictions calculated from tide height values which are outside of the 
range of values used for the initial regression analysis. For example, 
if all casts were acquired during a neap tide period, such as was the 
case for stations 2N and 2S, velocities predicted by the regression 
analysis for high or low tide heights of a spring tide period are subject 
to error. For stations 4SB and 4SA, predictions were calculated using 
a 9.0-foot tide height; station 4SB (Figure 3.2) was occupied during 
higher high tide, whereas station 4SA (Figure 3.3) was not occupied 
during this period. The second- and third-order analysis predicted 
velocities which probably were not physically possible at station 4SA. 
Therefore, a basic objective in using polynomial regression models for 
mean sound velocity predictions and corrections should be the 
determination of the maximum and minimum velocities expected, and the 
minimum range of tide where these velocities could occur. The most 
accurate method for identification of these values would be to measure 
the sound velocity variables at the same location covering more phases 
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of the lunar cycle. The danger of predicting outside of the observed 
range is illustrated by Figure 3.3, whereas Figure 3.2 shows that 
predicting within the observed range can be done accurately. 

Determination of the most efficient polynomial regression 
analysis to use for prediction of mean sound velocity was accomplished 
by examination of plots of the data arrays and comparison of the 
calculation of predicted velocities and confidence intervals provided by 
each analysis. 

(1) Data Array Plots. Examination of the data array 
plots indicated that the first-order (linear) analysis would be sufficient 
for most of the data arrays, although some of the arrays did suggest a 
curvilinear trend. Curvature could be detected in the plots of stations 
4NA, 4SB (Figure 3.2), 4SA (Figure 3.3), and 4NB (Figure 3.4), with 
maximum curvature in the plots of stations 4SA and 4SB. This may be 
an indication that the degree of curvature is proportional to the range 
of tide since data for stations 4SA and 4SB were acquired during the 
spring tide period, whereas data acquisition for stations 4NA and 4NB 
was carried out during a period of average tide range. This 
hypothesis cannot be verified without further investigation. 

(2) Predicted Velocities . A computer program was 

written to predict the mean velocity for five different heights of tide 
with the maximum and minimum heights being the same as those used in 
the regression analysis. The other three heights were equally spaced 
between the maximum and minimum. Rarely was there a significant 
difference (plus or minus 4 meters per second) between the predictions 
calculated by the three orders of polynomial regression (Table XIV). 
Mean velocities predicted using the the third-order regression analysis 
were in error for stations 2S and 5NB at the 13-meter depth. The 
cause of this error could not be explained. 

(3) Confidence I ntervals . Confidence intervals (Table 

XV) were computed for the predicted velocities using t values for 
90-percent confidence. In general terms, the confidence intervals were 
larger for all stations in the 1- and 3-meter depths, with the intervals 
becoming less with increasing depth. This corresponds with the lower 
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SOUND VELOaiY (M/5) 

1460.0 1465.0 1470 0 1475 0 1480 0 M85.0 H90.0 1495 0 1500.0 1505 0 



STATION 4SB 




O = 1ST ORDER PREDICTION CURVE 
X = 2ND ORDER FREDICTION CURVZ 
A = 3RD ORDER PREDICTION CURVE 

T 1 1 1 1 1 1 1 1 1 1 1 

- 2.0 - 1.0 0.0 1.0 2.0 3.0 4%0 5.0 6.0 7.0 0.0 9.0 10.0 

TDE ffiGHT (FT 



Figure 3.2 4SB--Prediction Curves--Depth = 5 m. 
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SOUND VELOQTY (M/5) 

HfeO.O M65.0 H70 0 K75 0 M80 0 1465,0 1490,0 K95 0 1500,0 1505 0 



STATION ASA 




-XO -1.0 0.0 ^0 XO 3.0 4sO B.O 6.0 7.0 9.0 0.0 10.0 

IDE l-BQ-iT (FT 



Figure 3.3 4SA--Prediction Curves--Depth = 5 m. 
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SOU^D VELOarY (M/5) 

0 M80.0 M85.0 VWO 0 M95 0 



STATION 4MB 











O 

18 

^-i 1 1 r* 

-10 -1.0 0.0 VO 



o = 1ST ORDER PFSacnON CURVE 
X = 2ND ORDER PREDICTION CURVE 
A = 3RD ORDER PREDICDON CURVE 

1 1 1 1 1 1 1 I I 

10 3.0 4%0 5.0 6.0 7.0 B.0 9.0 10.0 

TDE HQ6HT (FT) 



Figure 3.4 4NB--Prediction Curves--Depth = 5 m. 
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TABLE XIV 

Predicted Velocities--Depth = 5 m 



STN 

2N 


1.2 


TIDE 

2.6 


HEIGHT (ft) 

4.0 5.5 


6.9 


ORDER 

1 

2 

3 


1475.5 

1474.8 

1474.1 


1479.0 

1478.9 

1479.3 


1482.5 • 

1482.7 

1482.5 


1486.2 

1486.2 

1485.3 


1489 . 7 
1489.2 

1488.7 


STN 

2S 


0.3 


TIDE 

1.9 


HEIGHT (ft) 

3.6 5.2 


6.9 


ORDER 

1 

2 

3 


1476.8 
1479.5 
1479 . 3 


1480.0 

1479.6 

1479.8 


1483.3 
1481.9 
1481.9 • 


1486.5 

1486.1 

1486.0 


1489.9 
1492.8 
1493 . 1 


STN 

4NB 


0.0 


TIDE 

2.1 


HEIGHT (ft) 

4.1 6.2 


8.3 


ORDER 

1 

2 . 

3 


1472.8 

1474.4 

1474.4 


1476.6 

1475.-9 

1476.0 


1480.3 

1478.8 

1478.8 


1484.2 

1483.5 

1482.9 


1488.0 
1489 . 6 
1488.5 


STN 

4SB 


-1.5 


TIDE 

1.2 


HEIGHT (ft) 

3.8 6.5 


9.2 


ORDER 

1 

2 

3 


1465.0 

1469.2 

1469.5 


1471.5 
1471. 1 
1470.9 


1477.8 

1475.9 
1475.6 


1484.3 

1483.9 

1482.9 


1490.8 

1494.9 
1491.4 


STN 

5NB 


0.2 


TIDE 

2.2 


HEIGHT (ft) 

4.2 6.2 


8.3 


ORDER 

1 

2 

3 


1468.9 

1470.5 

1471.1 


1472.4 
1472.0 

1471.5 


1476.0 

1475.0 
1475.3 


1479.6 
1479.3 

1480.6 


1483 . 4 
1485.2 
1485.7 


STN 

5SA 


-1.6 


TIDE 

1.2 


HEIGHT (ft) 

4.0 6.7 


9.5 


ORDER 

1 

2 

3 


1460.2 
1461 . 6 
1461.9 


1463 . 1 
1462.7 
1462.5 


1465.9 

1465.0 

1464.7 


1468.6 

1468.3 

1467.2 


1471.5 

1473.1 

1468.8 


ORDER 


= Degree 


of polynomial regression 
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correlation coefficients found in the shallow depths. Larger confidence 
intervals were also found for data acquired at stations 2N, 2S, 4SA, 
and 4SB. These stations had lower correlation coefficients and higher 
residuals when compared to the other stations. For stations 2N and 2S, 
this may indicate less stability due to the close proximity of the open 
ocean. As stated before, data for stations 4SA and 4SB were acquired 
during spring tide conditions, which tends to indicate less stability for 
prediction. 

In most cases, increasing the degree of the 
regression analysis did not produce more accurate results in terms of 
confidence intervals. The standard error of estimate was sometimes 
lower with a higher order analysis, but the reduction in the degrees of 
freedom countered the effect of the lower standard error of estimate for 
computation of confidence intervals. 

Confidence interval results indicate that correcting 
soundings for mean velocity variation based on predictions provided by 
the regression analysis cannot be considered a complete success. 
Attaining the goal of predictions to within plus or minus 4 meters per 
second may not be possible in some cases, but the regression analysis 
does show a significant reduction in error when compared to only using 
average velocity and horizontal variation as the correction procedure. 

B. NOS BAR CHECK DATA 

Analysis of the bar check data revealed too many inconsistencies 
and required too much speculation for- comparison of the CREDDP data 
results. A major problem of the direct comparison method of the bar 
check is that examination of the sound velocity variables and any 
change in their magnitude is not possible. The bar check data did 
indicate the possibility of either an extremely high sound velocity 
variation or an instrument error which was not constant with depth. 

Most likely, the results show an undeterminable combination of sound 
velocity variation and instrument error. The problem of using the bar 
check data for any meaningful analysis is illustrated by Table XVI. 
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TABLE XV 

Confidence lntervals--Depth = 5 m 



STN 

2N 


1.2 


TIDE 

2.6 


ORDER 

1 

2 

3 


4.2 
4.1 

4.3 


3.9 

3.9 

4.0 


STN 

2S 


0.3 


TIDE 

1.9 


ORDER 

1 

2 

3 


4.3 

3.8 

3.9 


4.1 

3 . 6 

3.7 


STN 

4NB 


0.0 


TIDE 
2 . 1 


ORDER 

1 

. 2 

3 


3 . 1 
2.9 
3.3 


2.9 

2.7 

3.1 


STN 

4SB 


-1.5 


TIDE 

1.2 


ORDER 

1 

2 

3 


5.2 

4.9 

4.9 


5.0 

4.7 

4.7 


STN 

5NB 


0.2 


TIDE 
2 . 2 


ORDER 

1 

2 

3 


3.3 

3.0 

2.9 


3.2 

3.0 

2.8 


STN 

5SA 


-1.6 


TIDE 

1.2 


ORDER 

1 

2 

3 


3.9 

3.6 

4.5 


3 . 8 
3.5 
4.4 


units 


= ± m/s 





HEIGHT (FT) 



4.0 


IT) 

IT) 


6.9 


3.8 


3.8 


4.0 


3.7 


3 . 7 


3.9 


3.9 


3.9 


4.1 


HEIGHT (FT) 


3 . 6 


5.2 


6.9 


4.0 


4.1 


4.3 


3.5 


3.6 


3.7 


3 . 6 


3.7 


3.8 


HEIGHT (FT) 


4.1 


6.2 


8.3 


2.8 


2.8 


3.0 


2.6 


■ 2.7 


2.8 


3.0 


3.0 


3.2 


HEIGHT (FT) 


3.8 


6.5 


9.2 


4.9 


4.9 


5 . 1 


4.6 


4.6 


4.8 


4.6 


4.6 


4.8 


HEIGHT (FT) 


4.2 


6.2 


8.3 


• 3.2 


3.2 


3.3 


2 . 9 


2.9 


3 . 1 


2 . 8 


2.8 


2.9 


HEIGHT (FT) 


4.0 


6.7 


9.5 


3.7 


3.8 


3.9 


3 . 5 


3.5 


3 . 6 


4.4 


4.4 


4.6 
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TABLE XVI 

Bar Check Correction Data for Survey H-8422 



UPPER LOWER 



OBS 


MEAN 

ft 


LIMIT 

ft 


LIMIT 

ft 


RANGE 

ft 


DEPTH 

ft 


DEL V 
m/s 


88 


-1.0 


-0.6 


-1.4 


0.8 


5 


231.5 


90 


-0.8 


-0.2’ 


-1.4 


1.2 


10 


177.5 


90 


-0.7 


-0.1 


-1.4 


1.3 


15 


124.4 


90 


-0.7 


0.0 


-1.3 


1.3 


20 


93.5 


90 


-0.7 


0.1 


-1.4 


1.5 


25 


88.3 


87 


-0.7 


0.0 


-1.6 


1.6 


30 


78.4 


82 


-0.6 


0.0 


-1.3 


1.3 


35 


53.5 


73 


-0.5 


0.0 


-1.2 


1.2 


40 


44.4 


30 


-0.3 


0.6 


-1.1 


1.7 


45 


55.5 



DEL V represents the magnitude of velocity difference 
necessary to compensate for the range of the 
corrections if there were no instrument error. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Investigation of sound velocity variation in the Columbia River 
estuary revealed the following significant findings: 



Sound velocity variation in certain areas required the examination 
of how veloctiy changed with time and location if NOS accuracy 
standards were to be met. The procedure of computing only one 
corrector table based on an average velocity, or correction, was 
not acceptable. 



CREDDP data analysis showed a strong relationship between tide 
height and sound velocity, which stimulated the development of a 
method to account for the variation of sound velocity with time. 
Using a regression analysis, the accuracy of velocity correctors 
was improved significantly, although there were still instances 
where NOS standards could not be met. 



NOS bar check data analysis showed a significant variation in 
corrections, although the variation did not correlate with tide 
height as the CREDDP data did. Analysis of the variation found 
at each 5-foot increment of depth revealed inconsistencies which 
could not be explained. 



B. RECOMMENDATIONS 



The following recommendations are made: 



Further study of sound velocity variation must be carried out for 
verification of ideas and methods presented in this thesis. 

Analysis of data acquired over longer time periods at the same 
locations will allow a more complete study of the variation 
associated with different lunar phases and the effect of changes in 
river flow. NOS CTD data were acquired in the Columbia River 
estuary during the same time period as CREDDP; these data may 
allow study ot the sound velocity variation over longer time 
periods . 



NOS bar check data analysis showed that the hydrographer shpuld 
be equipped with an efficient method of measuring sound velocity 
in addition to determining instrument error. Knowing the value of 
sound velocity, detection of significant changes and trends is more 
efficient and accurate. Bar checks should not be used as the 
primary method of correcting soundings, but should be reserved 
for determination of instrument error. 



51 



Hydrog rap hers conducting surveys in areas where reversing tides 
could cause a significant variation in sound velocity may benefit by 
noting the height of tide for each cast or bar check acquired for 
determining sound velocity corrections. If significant variation is 
detected, the possibilty of predicting and correcting for temporal 
variation using the tide height as an independent variable should 
be investigated. FORTRAhT computer programs written for the 
methods presented in this thesis are listed in Appendices A, B, C, 
D, E, and F. 



In areas where temporal variation of sound velocity is significant, 
the hydrographer may benefit by using a sound velocimeter to 
correct soundings. The velocimeter would require only a few 
minutes to deploy and retrieve. The hydrographer could determine 
sound velocity periodically throughout tne working day without 
significantly hindering operations. 
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FORTRAN PROGRAM PACKAGE- - SOUND VELOCITY COMPUTATIONS 













-'c 










•Jc 










-Jc 












LO 










Z 








-Jc 


O 






•)% 


-)% 


H 






*■){ 




H 










<J 










H 










P 








-j? 


PM 


w 








S 


Q 








O 


P 


00 




-jc 


CJ 


H 


GN 


•){ 








rM 












-j? 




H 


• 


X 


“)C 


-)< 


M 


Q 


CJ 




•j? 


U 




piJ 


-J' 


•)% 


o 


3: 






-jc 




P 


5; 


•){ 


•)< 


w 


o 






•>? 


> 


p 


















Q 










2: 






-)C 




D 










O 










LO 








•j: 


















.’j 










-;c 














*:? 









PI5 




3: 








5/1 


• 


s 


pp 


o 


o 


M 


p 


H 


H 


<J 


W 


D 


P 


C/ 




W 


P 




C 


c/3 


P 


- 


M 


JS 


O 


O 


M 


c/3 


W 


X 


PM 


M 


<; 




H 


o 




2: 


CxJ 


M 


p 


c/3 


H 


D 


O 




P 


H 


M 


M 


C/3 


CJ 


P 


O 




X 




MJ 


H 


> 


M 




CJ 


Q 


O 


2: 


P 


D 


u 


O 


> 


c/3 


Q 


C/3 


P 


W 


P 


H 


O 


P 


c/3 


Pm 




S 


CJ 


O 


M 


u 


P 



o 




o 

PM O 



M < 



UCJCJUCJUUCJ 



u u 



CJCJCJCJCJUUUUCJ 



53 



IMPLICIT REAL-4 (A-H,0-Z) 
CUARACTER-'51 H(8) 
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READ(4,70) H(J) 

FORMAT (A51) 

CONTINUE 

READ(4,80) H1,DATE,H2,TH,H3 ,TM 
FORMAT ( IX, A13 , All , 3X , A12 , IX , 12 , A1 , 12 ) 



THO=TH+0.0 

TMI=( (TM+0.0)/10.0)/6 .0 

TIME=THO+TMI 

WRITE (7, 75) 

FORMAT ( / / , IX , ’ STATION ’ , 7X , ' DATE ’ , 8X , ' TIME ' 
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VP(J)=P(J)*( . 160272+P(J)*(1.0268E-5+P(J)*(3.5216E-9-P(J)*(3.36 
03E-12)))) 

VS(J)=(S(J)-35 . )*(1. 39799+ (S(J) -35. )*1 . 69202E- 3 ) 

A1(J)=T(J)*(- .011244+3 .1580E-8*P(J) +T(J)*( 7 .7711E- 7+1. 5790E-9" 

P(J))) 



A2(J)=P(J)*(7.7016E-5-1.2943E-7*P(J)) 

A3(J)=P(J)*T(J)*(-1.8607E-4+T(J)*(7.4812E-6+T(J)*(4.5283E-8))) 

A4(J)=P(J)*P(J)*(T(J)*(-2.5294E-7+T(J)*1.8563E-9)+P(J)*T(J)*(- 

1.9646E-10)) 
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READ(4,110) DT(I) 

IF(I.EQ.l) GO TO 25 

IF(DT(I) .LT.DT(I-l)) DT ( I ) =DT ( I ) + 24 . 
DT(I)=DT(I)+24.0 
110 FORMAT ( IX, F5. 2) 



BEGIN LOOP FOR COMPUTATIONS 
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IMPLICIT REAL-4 (A-H,0-Z) 
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110 F0RMAT(1X,A3,2X,F5. 1,2X,F5. 1,2X,F7. 1,2X,F5. 1,2X,F5.2) 

5 CONTINUE 

STOP 
END 
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FORTRAN PROGRAM PACKAGE- -FIRST-ORDER REGRESSION. 
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INTEGER B,I,M 
CHARACTER*! STN(99) 

DIMENSION HT(50) ,V(50) ,ER(50) 
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= B+1 

=C+HT(I) 

=D+ (HT(I)**2) 
=E+V(I) 



Q 

c 



PQ U 



Q 






O 

M 

H 

C 

D 

O* 

CxJ 

S 

o 



CO 

•J 

< 

D 

Q 

H 

CO 

pit 

O 




CM 

-)c 



M 

H 

CJ 

H 

Q 

W 

PU 

PX4 

o 



CM 



s 

CO 



CO 

hJ 

C 

D 

Q 



C 

Q 











a 




CO 




00 






CO 






C 






M 




V 












'W' 




H 




PX4 






Pl] 






1 






CO 




CM 










• 






2: 


<d 








Pi 












PQ 




• 










CM 


1 




Px:i 




X 


















Pi 


/"“V 


\D 










C 


PQ 




M 


CM 


m 






P4 






M 


/— N 


CM 




M 


U* 
















U 


<d 








O 








M 




Q 


V— ✓ 


r* 










Q 


— 




M 




CM 












H 


'W' 


•/* 


X 


Pi 


XI 




/--~s 






X 


(~i 






/-“N 


• 






X 






X 


Pi 




< 


PQ 


r* H 










< 










un 






o 


• 




>* 

y* 


CJ 


M 






r-l 




M 












W 




P&4 






M 


o 




CM 


N«-^ 


N— ✓ 


Pi 


M 


•> 




^ " 






r— 1 






o 


X 


#« 






H 


PJ 




< 


CO 


Pi 


CJ 




CM 




> 






C 


/<— N 


Q 


u 


H 


X 






C 


Pi 




"W ^ 


PQ 


PQ 


PQ 


/— s 


M 










Q 






CO 


on 








c 




1 


C 




S 


O 












CO 


— 




Q 










D 


D 




^ — V 


+ 


+ 




r—l 






H 






H 






S 


/■ — s 


cn 






CJ 


CJ' 




M 


r—l 


CM 


s 


r— ^ 


O 










Z 


O 


CJ 


c 


un 


.<< 








CO 


PQ 




X 


S 






r—l 




H 






W 






CM 








c 






> 


D 


D 


s 


r*^ 


^ ^ 




X 






M 


U 


r-^ 


Q 


r-^ 


>< 






u 


CO 


f—*i 


II 


CO 


CO 


o 


'W' 


H 






UJ 


PjJ 


CJ 




< 


M 






O 


o 






II 


/■ V 


tl 


II 


M 


PQ 






+ 


:d 


Z) 


M 


>*-✓ 








V ^ 


• 


• 


Pi 


C 


M 


M 




CM 


H 


H 


y 




O 


s 


3; 




1 


1 


X 




H 


o 


o 


o 


D 






s 


S 


C 


H 


2 




It 


H 


M 




P£] 


CJ 


o 






n 


II 




Q 


O 


Pi 




J) 


> 


Pi 


o 


PQ 


o 


H 


H 


PxJ 


'W 


n_^ 


M 


H 






CM 




M 


CO 


PQ 


CO 


CO 


Pi 


3: 


Cx. 


J) 




s: 


:s 


O 


M 


tl 


H 


M 


2 


s 




P4 


CO 










PQ 






s 




o 


o 


u 


r— 1 


CM 


C 


Pi 


o 


:j 




o 


PQ 


O 








CO 






M 




u 


u 




c 


< 


H 


::^ 


pu 


CO 


CO 


o 


Pi 


PQ 








PQ 






H 








UJ 






CO 






















O 






Z 








H 


















p4 
















O 








D 






H 










s 


O 










H 






CJ 








P-» 






IS 










M 












S 














S 






M 










CJ 


s 










M 














o 






Pi 










PJ 


D 










Pi 














CJ 






P-» 




un 






PQ 


CO 










P4 




o 






o 


o 












CM 
























o 






CM 


CJ 






CJ 




r—l 






CJ 


CJ 










CJ 






CO 



70 



PRINT SUM OF RESIDUALS, SUM OF RESIDUALS SQUARED, AND NUMBER OF OBS . 
WRITE (7, 125) SUM1,SUM2,B 
125 FORMAT(/ ,28X,F7.2,3X,F7.2,3X,I2) 

I 

5 CONTINUE 
STOP 
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FORTRAN PROGRAM PACKAGE- -SECOND-ORDER REGRESSION. 
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VARIATION OF SOUND VELOCITY. 



INTEGER B,K 
CHARACTER*3 STN(50) 

DIMENSION HT(50) ,V(50) ,ER(50) 
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READ IN TIDE HEIGHTS AND OBSERVED SOUND VELOCITIES. A FLAG OF 



999,99 IS TO BE PUT AFTER LAST VELOCITY TO SIGNAL END OF LOOP. 
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A2=(G-((F*J)/H)>(A1*(((D^^F)/H)-C)))/L 
A3=(J- (A1*D)- (A2*F) )/H 

PRINT STATION ID AND COEFFICIENTS OF PREDICTION EQUATION. 
V/RITE(7 ,120) STN,AI,A2,A3 
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FORTRAN PROGRAM PACKAGE- -PREDICTED VELOCITY/ CONFIDENCE INTERVAL 
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REAL V(99) ,T(99) ,S(99) ,H(99) ,HM(99) ,SH(99) ,SV(99) ,SH2(99) ,SV2(99) , 
*A(99) ,B(99) ,LC(99) ,UC(99) ,T1(99) ,WI(99) ,S2(99),CI(99),A2(99) ,B2(99 
*) ,C2(99) ,T2(99) ,V2(99) ,LC2(99) ,UC2(99) ,CI2(99) ,SUM2(99) ,SUM(99) , 
“DEP 

INTEGER N(99) ,I ,K, J,D,K2,L 

































































































PJ 






































> 






































H 














• 
























H 














CO 
























X 














X 
























o 






p4 








c 
























PJ 






O 








X 
























CO 














PQ 
























X 






CO 








M 
























o 






PD 








CO 
























o 






PP 








P 






























C 








P 
























PU 






X 
































PD 






O' 








P 
























M 






CO 








o 












































































CO 






























o 








p 
























X 














PP 






























S 








C 


















• 












X 








X 


















CO 






PD 






CO 






CSJ 


O' 




N 














'Z 






PQ 










M 


M 


CO 




CSJ 














o 






















/— s 


• 














H 






O 






CO 




X 


>P 


P 


M 


vO 














H 






H 






H 






CSJ. 


o 


'« — ' 


p 


























X 










CSJ 
















> 






X 






PD 


• 


M 


CSJ 


X 


X 


X 














pp; 






O 






M 


CO 


'w' 


• 


X 


X 


CSJ 




















M 






U 


X 


X 


vO 


CO 


CO 
















CO 






H 






M 


o 


X 


Pi 






CSJ 














P3 






C 






p4 


M 


CO 




PQ 




• 














O 






H 






p4 


H 




X 


X 


M 


to 




















CO 






PD 


C 




csl 


c 


>«—✓ 


p 


























O 


> 


M 






CSJ 
















u 






PD 






U 


P^ 


N <' 


csj 


CO 


O 


X 














CxJ 






•X 








PD 


PQ 


• 


H 


- 


CSJ 




• 










X 






o 






P^ 


CO 




LO 


X 


/— s 






CO 










u 












PD 


PQ 




p 


P 


M 


CSJ 




w 
















X 






PP 


O 


M 




H 




• 




D 










ppi 






:< 






pi; 






>p 


O 


CSJ 


LO 














c 






X 






o 


P4 


c 


CSJ 


M 


X 


p 




C 










PQ 






H 






1 


o 






P 








> 






















H 




N 


CSJ 


P 




X 








H 






ppi 






PP 






CO 


pp 


M 


• 


P 


M 


CSJ 








"■ — ^ 






O 






pp; 


X 




pi; 


PD 


'' — ' 




O 


"■ — ' 




ON 


















o 


o 




M 


X 


H 


P 


o 


CSJ 


CSJ 


. 


PQ 




H 






H 






X 


M 




p^ 


s 


CO 






C 


• 




C 










CO 


pu 






H 






X 




x; 


pp 






H 


H 




N 


m 




C 


C£J 




CO 


C 






X 




m 


p 


s 


P 


CO 






o 


• 




U 


o 




::5 


H 




pp 




o 




X 


1—1 






H 


o 


c\J 


uO 










o 


CO 




M 


PQ 


o 


m 


pp 


o 


X 


m 


1 


m 








Pz-i 




r-l 


rJ 




r-l 




X 


I— 1 


C 


o 


1— 1 


vO 




H 






'w' 




O 


CN 


• 


rJ 


CO 


T-^ 


X 






Vw' 


1 








2S 




m 


H 






CSJ 


<r 


<< 


PD 




o 




r-l 


H 


PQ 


1— 1 


H 


u 


UJ 


n 




<] 




X 


r— 1 


Pm 




H 


r-l 


H 








X 






H 


Q 


H 


Q 




Z) 


H 




■>w/ 


Pm 


<J 


II 


H 


CO 


PQ 




o 


X 




U 


ZD 




C 


p2 


2= 


Pu 


<r 


H 


o 


X 


M 




X 


< 


2 


o 


< 


2 


< 


H 


O 




o 


M 


UD 


'W' 




o 


o 




H 


< 


PD 


o 


p 


P 


o 




CO 


r-l 


Pil 


u* 


H 


O 


pp 




rJ 


M 


LO 


CO 


X 


PP 


p 


CO- 


pp 


p 












S 




<< 






CO 






PQ 














o 


o 






O 


O 


C£1 


O 


CO 


PD 


O 


p:i 


M 






Q 






U 




o 






O 


< 


PP 


PL4 


H 


Q 


Q 




CO 






< 


















CjJ 






X 






PD 


PD 






P 


















p:: 






H 


M 




P^ 


PP 






pp 














o 








OS 














o 






T— 1 










csj 


o 






CSJ 














o 






o 










r-l 


1—1 






r-l 














p 






1— 1 




o 










u 






U 


a 




o 


U 






u 







77 



READ AVERAGE TIDE HEIGHT, STANDARD DEVIATION OF TIDE HEIGHT, AND 



STANDARD DEVIATION OF OBSERVED VELOCITIES. 

READ(2,110) HM(I) ,SH(I),SV(I) 

110 FORMAT ( 14X , P4 . 2 , 15X , F4 . 2 , 4X , F4 . 2 ) 

READ IN TIDE HEIGHTS FOR WHICH PREDICTIONS ARE TO BE BASED. 
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PRINT STATION ID AND TIDE HEIGHTS. 



211 F0RMAT(//,1X, 'STN' ,15X, 'TIDE HEIGHT (FT ) ’ , / , IX , A3 , 3 ( 9X , F6 . 1 ) , / ) 

BEGIN LOOP FOR PREDICTED VELOCITIES AND CONFIDENCE INTERVALS. 
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